A Critical Assessment of Nonlinear Force-Free Field Modeling of the Solar Corona 

for Active Region 10953 
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ABSTRACT 

Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics 
and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that 
NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically 
constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models 
have been found to have markedly different field line configurations and to provide widely varying estimates of the 
magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force- 
free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is 
dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several 
major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active 
Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the 
process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic 
field data covering larger areas are needed so that more electric currents associated with the full active regions of interest 
are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, 
and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better 
transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of 
the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems. 
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1. Introduction 

The structure and evolution of the magnetic field (and the 
associated electric currents) that permeates the solar atmo- 
sphere play key roles in a variety of dynamical processes 
observed to occur on the Sun. Such processes range from 
the appearance of extreme ultraviolet (EUV) and X-ray bright 
points, to brightenings associated with nanoflare events, to the 
confinement and redistribution of coronal loop plasma, to re- 
connection events, to X-ray flares, to the onset and liftoff of 
the largest mass ejections. It is believed that many of these ob- 
served phenomena take on different morphologies depending 
on the configurations of the magnetic field, and thus knowl- 
edge of such field configurations is becoming an increasingly 
important factor in discriminating between different classes 
of events. The coronal topology is thought to be a critical 
factor in determining, for example, why some active regions 
flare, why others do not, how filaments form, and many other 
topics of interest. 

One model of the coronal magnetic field B assumes that 
the corona is static and free of Lorentz forces, such that J X 
S = 0, where J = c V X S/47r is the current density. This 
means that V X S = olB, and thus any electric currents must 
be aligned with the magnetic field. Because V • ^ = 0, it can 
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be shown that S • Va = 0, demonstrating that a is invariant 
along field lines of B. The scalar a is in general a function 
of space and identifies how much current flows along each 
field line. In cases where a varies spatially, the problem of 
solving for B (and ol) is nonlinear. Solving for such nonlinear 
force-free fields (NLFFFs) requires knowledge of B over the 
complete bounding surface S enclosing the solution domain. 
To be compatible with a force-free field, it is necessary for 
these boundary data B\s satisfy a number of consistency 
criteria, which we outline in ^ an d which are expl ained in 
detail in lMolodenskiil (Il969h and in lAivl(ll984ll989k . 

In analyzing solar active regions, localized maps of the 
photospheric vector field are typically used for the lower 
bounding surface B\z^, and potential fields are used for the 
other surfaces. (For the Cartesian models discussed herein, 
we use the convention that the z axis is normal to the pho- 
tosphere, which is located at height z = zq.) The availabil- 
ity of vector field maps produced by recent instrument suites 
such as the Synoptic Optical Long-term Investigations of the 
Sun (SOLIS) facility and the Hinode spacecraft, building on 
earlier work done in Hawai'i with data from the Haleakala 
Stokes Polarimeter (HSP) and by the Imaging Vector Mag- 
netograph (IVM) as well as from the HAO/NSO Advanced 
Stokes Polarimeter (ASP) at Sacramento Peak in New Mex- 
ico, has spurred investigations that employ coronal-field mod- 
els based on such measurements. We anticipate that such re- 
search will intensify when regular, space-based vector field 
maps from the Helioseismic and Magnetic Imager (HMI) in- 
strument on board the Solar Dynamics Observatory (SDO) 
become available. 

One goal of NLFFF modeling is to provide useful esti- 
mates of physical quantities of interest (e.g., connectivities, 
free energies, and magnetic helicities) for ensembles of ac- 
tive regions, so that these active regions may be systemati- 
cally analyzed and intercompared. The use of static, force- 
free models mitigates some of the computational difficulties 
associated with solving the more physically realistic, time- 
dependent problem, as running such dynamical models at the 
desired spatial and temporal resolutions for multiple active re- 
gions typically exceeds current computing capabilities. 

There exist several previous studies of individual active re- 
gions where NLFFF models are shown t o be compatible with 

various stru ctures in t he corona (e.g., Regn ieret"aD 2002 ; 

Regnier & A mari 20 04l:lwiegelmann et al.l2 005: Reg nier & Canfieldl 
20061: ISchriiver et akboO SV Several of these studies provide 
evidence of good alignment between NLFFF model field lines 
and the locations of observed features such as coronal loop 
structures observed in EUV and X-ray images. Others show 
that the locations of sigmoids, twisted flux ropes, and/or field 
line dip locations coincide with analogous features in the 
NLFFF models. Such studies are certainly encouraging, but 
still it remains difficult to conclusively determine whether 
these models match a significant fraction of the coronal mag- 
netic field located within the volume overlying an entire active 
region. 

As part of a long-lasting (e.g.. lSakuraill98ll:lMcClymont et al. 



1997|) effort to develop methods that generate more robust 



NLFFF models, a working group (in which all of the au- 
thors of this article are participating) has held regular work- 
shops over the past several years. The previo us resu lts from 
this collaboration a re pr esented in ^S chriive r et al.l (I2OO6), 
Metcalf etall (l2008k . and ISchrijver et al.1 (120081) . Since the 
launch of Hinode in 2006, we have applied multiple NLFFF 
modeling codes to a few active regions for which Hinode 
vector magnetogram data are ava ilable and for which n onpo- 
tential features are evident (e.g.. ISchrijver et aDl2008 ). The 
resulting NLFFF models generally differ from each other in 
many aspects, such as the locations and magnitudes of cur- 
rents, as well as measurements of magnetic energy in the so- 
lution domain. In this article, we identify several problematic 
issues that plague the NLFFF-modeling endeavor, and use a 
recent Hinode case to illustrate these difficulties. We describe 
one representative data-preparation scheme in ^ followed 
in by a comparison of field lines in the resulting NLFFF 
models to two- and three-dimensional coronal loop paths, the 
latter determined by analyzing pairs of stereoscopic images. 
In ^ we explain the primary issues that we believe to impact 
our ability to reconstruct the coronal field in a robust manner, 
and also identify and discuss the alternate data-preparation 
scenarios we tried in addition to those presented in Con- 
cluding remarks are presented in ^ 

2. Construction of NLFFF Models for AR 10953 

Several NLFFF extrapolation algorithms (each implement- 
ing one of the three general classes of extrapolation methods) 
were applied to boundary conditions deduced from a scan 
of NOAA Active Region (AR) 10953, taken by the Spectro- 
Polarin ieter (SP) in strument of the Solar Optical Telescope 
(SOT) (ITsuneta et al. 2008) on board the Hinode spacecraft. 
The Hinode/SOT-SP scan of this active region started at 
22:30 UT on 2007 April 30 and took about 30 min to com- 
plete. As the scan progressed, polarization spectra of two 
magnetically sensitive Fe I lines at 630L5A and 6302. 5A 
were obtained within the 0'.'16x 164'' slit, from which Stokes 
IQUV spectral images were generated. For this scan (in "fast- 
map" mode), the along- slit and slit- scan sampling was 0'.'32, 
and the total width of of the scan was 160''. AR 10953 pro- 
duced a C8.5 flare about two days after this Hinode/SOT-SP 
scan, and a C4.2 flare about four and a half days after this 
scan, but otherwise the active region was flare-quiet above 
the Cl.O level. Jmages from the X-Ray Telescope (XRT) 
on board Hinode around this time show 



(iGolub etal.ll2007h 
a series of bright loops in the central region of AR 10953 

(Fig. [a). 

The NLFFF algorithms need vector magnetic data as 
boundary conditions, and determining these boundary maps 
comprises the first step in constructing a NLFFF model. The 
conditions pertaining to the lower boundary are determined 
from a map of the photospheric vector magnetic field from 
the Hinode/SOT-SP instrument. The magnetic components 
parallel to and transverse to the line of sight, 5los and Bi, 
are functions of the circular and linear polarization signals, 
respectively. Constructing B\zq requires assuming an at- 
mospheric model (in this case Milne-Eddington) and de- 




Fig. 1. — A series of coaligned images of AR 10953 (with the same 10° gridHnes drawn on all images for reference), (a) Time- 
averaged and logarithmically scaled Hinode/XRT soft X-ray image, and (b) with the best-fit Wh~ model field lines overlaid, (c) 
STEREO-A/SECCHI-EUVI 171 A image, (d) Trajectories of loops, as viewed from the perspective of an observer located along the 
Sun-Earth line of sight and determined stereoscopically from contemporaneous pairs of images from the two STEREO spacecraft, (e) 
Same visualization as panel (d) but viewed from the side. The solid black cubes in panels (d) and (e) outline the full 320 x 320 x 256- 
pixel NLFFF computational domain, and the interior dotted black square outlines the base of the smaller 1 60 x 160 x 160-pixel 
volume (covering most of the Hinode/SOT-SP scan area) used for the field line maps of Fig. [2] and for the metrics in Table [T] The 
STEREO-loop points are colored blue if outside the NLFFF computational domain, or are colored according to their misalignment 
angle (j) made with the field lines from the Wh~ solution. Yellow is indicative of < 5°, red of ^ > 45°, with a continuous 
progression from yellow through orange to red for 5° < (j) < 45°. On the bottom face of the large cube is displayed the map 
used during the NLFFF modeling, which includes higher-resolution data from Hinode/SOT-SP embedded in SOHO/MDI full-disk 
magnetogram data. The magnetogram images saturate at ±1500 Mx cm~^. 



termining which combinations of magnetic field strengths 
and filling factors produce th e observed polarizatio n sig- 
nals (e.^.. ISkumanich & LitesI [198 7: Klimchuk et J] ll992|; 
iBorrero et al.ll2007l) . 5los has uncertainties that are typically 
an order of magnitude less than Bt. 

The next step involves removing the ambiguities in the 
components of Bt that arise due to the property that the same 
linear polarization signal can be produced by either of two 
magnetic field vectors differing by 180° of azimuth in the 
transverse plane. We choose to perform the disambiguation 
using the interactive Azimuthal Ambiguity Method (AZAM), 
which is one of several metho ds have been devised and tested 
to resolve this ambiguity (see iMetcalf et al.|[2006l and refer- 
ences therein). 

After disambiguation, the B\zq map for AR 10953 is used 
to produce potential field data with which the extrapolation 



codes will initialize the computational domain. Our approach 
is to specify the computational domain (having an enclosing 
surface S) that contains much of the coronal volume overly- 
ing the active region of interest, such that the lower bound- 
ary includes the area for which vector magnetogram data are 
available. The initialization field is calculated by embedding 
the Hinode/SOT-SP vector magnetogram data in a larger line- 
of-sight magnetogram obse rved by the Michelson Doppler 
Imager (MDI) instrument ( Scherrer et al.l ll995) on board the 
Solar and Heliospheric Observatory (SOHO) spacecraft (as 
shown in Fig. [TJi). Then, the potential field coefficients cor- 
responding to this enlarged footprint are determined, from 
which the potential field in the 320 x 320 x 256-pixel NLFFF 
computational domain is computed. In addition, the vector 
field boundary conditions for the side and top boundaries of 
the computational domain are taken from this same poten- 



tial field extrapolation, primarily because we expect that the 
coronal magnetic field becomes largely potential away from 
the center of the active region, but also because it is useful 
to specify how unbalanced flux emanating from this active 
region connects to flux of the opposite polarity located else- 
where on the Sun. 

The embedded lower-boundary data are then sampled onto 
a uniform, helioplanar, 320 x 320-pixel grid having 580 km 
pixels, such that the footprint of the computational domain 
spans a 185.6-Mm-square area. The region for which Hinode 
vector magnetogram data for AR 10953 were available com- 
prise about a lOO-Mm-by-115-Mm subarea of the full lower 
boundary footprint, outside of which the horizontal compo- 
nents of B\zQ are set to zero. Thus, in this peripheral re- 
gion outside the Hinode/SOT-SP field of view, the field on 
the lower boundary can either be considered as purely vertical 
(for force-free methods which use all three components of the 
field as boundary conditions), or equivalently as having zero 
vertical current density (for methods which use the vertical 
component of the field together with the vertical component 
of the current density). 

Next, to be consistent with a force-free field, it is neces- 
sary (but not sufficient) that the entire boundary fie l d B\q 
satisfy seve ral cri t eria, a s delineated in iMolodenskiil (|l969|) 
and in Alv (11984 [19891) : namely, (1) the volume-integrated 
Lorentz force must vanish, (2) the volume-integrated mag- 
netic torque must vanish, and (3) the amount of negative- 
polarity flux through having a given value of a must 
equal the positive-polarity flux through with this same 
value of a. The first two criteria are relations involving var- 
ious components of ^Is-, and are derived from volume in- 
tegrals of the Lorentz force and its first moment. The third 
("^-correspondence") relation operates over all values of a 
present on ^l^. 

There is of course no guarantee, however, that the values 
of B\zQ, coupled with the potential field of B for the comple- 
ment of the enclosing surface, together satisfy these consis- 
tency criteria. Our working group attempts to deal with this 
problem by preprocessing the boundary data before feeding 
them to the extrapolation codes. The preprocessin g scheme 
used here (developed by Iwiegelmann et al.l |2006|) seeks to 
adjust the components of B\zq so as to satisfy the first two 
consistency criteria while minimizing the deviations of B\zq 
from their measured values. During this preprocessing step, 
spatial smoothing is also applied to B\zq to attenuate some of 
the small-scale magnetic fluctuations that likely die off shortly 
above the photosphere. 

Finally, we apply the various NLFFF algorithms to these 
boundary and initial data. Several methods for calculating 
NLFFF models of the coronal magnetic field have been de- 
veloped and implemented in recent years, including (1) the 
optimization method, in which the solution field is evolved 
to minimize a volume integral such that, if it becomes zero , 
the field is divergen ce- and force-free (Wh eatland et aLlEoOOl: 
Wiegelmannl 12004 ) : (2) the evolutionary magnetofrictional 
method, which solves the magnetic induction equation us- 
ing a velocity field that advances the solution to a more 



force-free state dVang et al.lll"986[ IValori et al.ll2007k : and (3) 
Grad-Rubin- style current-field iteration procedures, in which 
currents are added to the domain an d the magnetic field 
is recomputed in an iterative fashion ( Grad & RubinI [l958[ 
Amari et al.ll2006l:lwheatlandll2006h . Some of these methods 
have been implemented by multiple authors. For brevity, we 
omit detailed explanations of these numerical schemes as im 



plemented here and instead dire ct the reader to lSchrijver et al 



12006) andlMetcalf et al. I (l2008h . and references therein. 

Although th ese methods work we ll when applied to sim- 
ple test cases dSchrijver et al.l l2006l) . we have found that 
the results from each of the methods typically are not con- 
sistent with each other when applied to solar data. The 
resulting magnetic field configurations differ both qualita- 
tively (e.g., in their connectivity) as well as quantitatively 
(e.g., in the amount of magnetic energy contained within 
th em). In discussing th e results from the solar-like test case 
of iMetcalf et al.l (|2008|) . we described some likely causes of 
such discrepancies amongst the models. In what follows, we 
illustrate these problems in greater detail using the (solar) data 
set at hand. 

3. Comparison with XRT and STEREO Loops 

The results of twelve extrapolations for AR 10953 (includ- 
ing the potential field), based on the data-preparation steps 
described in are summarized in Table [T] and Figure [2l 
Table [T] contains domain- averaged metrics characterizing the 
center of the active region (corresponding to the region sur- 
rounding the leading, negative-polarity sunspot), and Figure[2] 
shows representative field lines in this same volume for each 
of these models. This central region is a 160 x 160 x 160-pixel 
volume, chosen to cover the portion of the lower boundary 
containing much of Hinode/SOT-SP magnetogram data (i.e., 
where we have some knowledge about the currents passing 
through the photosphere), and is fully contained within the 
larger 320 x 320 x 256-pixel computational domain. 

The models considered in Table [T] and Figure [2] are the 
current-field iteration method as run by Wheatland using the 
values of a in either the negative or positive polarity (here- 
after "Wh~" and "Wh+", respectively) ; the finite-elemen t 



Grad-Rubin-style method (FEMQ in lAmari et al.l l2006h 



run using two different parameter sets by Amari ("Ami " 
and 'Am2~"); the vector-potenti al Grad-Rubin-like method 
(XTRAPOL in lAmarietal1l200d) by Canou ("Can""), or by 
Regnier using the values of a in either the positive ("Reg+") 
or negative ("Reg~") polarity; the optimization method using 
grid refinement as run by Wiegelmann ("Wie") or McTiernan 
("McT"), or no grid refinement as run by Thalmann ("Tha"); 
the magnetofrictional method using grid refinement as run by 
Valori ("Val"); and the initial potential solution ("Pot"). 

We find that the Ami", Am2~, Can", and Wh~ current- 
field iteration models contain between 18% and 25% more 
energy than the potential solution, and have smaller residual 
Lorentz forces and smaller average V • B than the other mod- 
els. In addition, the Ami", Am2~, and Can" models find a 
strongly twisted flux rope in equilibrium, whose foot points 
are anchored southeast of the main spot (mostly outside of 



Table 1 

NLFFF Model Extrapolation Metrics^ for AR 10953 
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(CW sm 6) " 


(xlO^) 


1 M\ f 


Pot 


1.00 




0.02 


24° 


Wh+ 


1.03 


0.24 


7.4 


24° 


Tha 


1.04 


0.52 


34. 


25° 


Wh- 


1.18 


0.16 


1.9 


27° 


Val 


1.04 


0.26 


71. 


28° 


Aml- 


1.25 


0.09 


0.72 


28° 


Am2- 


1.22 


0.12 


1.7 


28° 


Can- 


1.24 


0.09 


1.6 


28° 


Wie 


1.08 


0.46 


20. 


32° 


McT 


1.15 


0.37 


15. 


38° 


Reg- 


1.04^ 


0.37 


6.2 


42° 


Reg+ 


0.87^ 


0.42 


6.4 


44° 



^AU metrics were evaluated over a 1 60 x 160 x 160-pixel comparison vol- 
ume (whose base overlaps much of the Hinode/SOT-SP scan area and is shown 
as a dotted line in Figs. [TJi,e), with the exception of (</)}, for which the full 
320 X 320 X 256-pixel computational domain was used. The models are listed in 
order of (cf)). 

^As listed in the models are the initial potential solution ("Pot"); the 
current-field iteration method as run by Wheatland using the values of a in 
the negative ("Wh~") or positive ("Wh+") polarity; the finite-element Grad- 
Rubin-style method as run by Amari ("Ami"" and 'Am2~"); the vector- 
potential Grad-Rubin-like method by Canou ("Can~"), or by Regnier using 
the values of a in the negative ("Reg~") or positive ("Reg+") polarity; the 
optimization method using grid refinement as run by Wiegelmann ("Wie") or 
McTiernan ("McT"), or no grid refinement as run by Thalmann ("Tha"); and 
the magnetofrictional method using grid refinement as run by Valori ("Val"). 

'^E/Epot is the total magnetic energy relative to the initial potential field 
solution for the comparison volume. 

^The (CWsin^) metric is the current-weighted average of sin^, where 6 
is the angle between B and J in each model (with 0° < < 180°). For 
perfectly force-free fields, (CWsin^) = 0. 

^The {\fi\) metric is the mean over all pixels i in the comparison volume of 
the absolute fractional flux ratio \fi\ = |(V • B)i\/(6\B\i/ Ax), where Ax 
is the grid spacing. The {\fi\) metric is a measure of how w ell V • S = is 
satisfied in the models (cf. eq. [15] o f lWheatland erS]l2000h . with divergence- 
free fields having {\fi\) = 0. 

■^The quantity (0} is the mean difference in angle between the stereoscopi- 
cally determined loops and the NLFFF model field lines (with 0° < < 90°), 
averaged over the full NLFFF computational domain. 

^The Reg~ and Reg+ solutions use closed boundary conditions for the side 
and top surfaces through which no magnetic flux is transmitted, and thus are as- 
sociated with a different potential field than the Pot solution. When comparing 
the Reg~ and Reg+ solutions to the potential field associated with these closed 
boundary conditions, the values of E/Epot are 1.23 and 1.04, respectively. 
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Fig. 2. — Representative field fines in the central portion of the active region for each NLFFF model listed in Table [TJ The cubes 
shown here comprise the same 160 x 160 x 160-pixel subvolumes excerpted from the full 320 x 320 x 256-pixel computational domain. 
(The base of this sub volume is the region indicated by dotted lines in Figs.[T]l,e.) The starting locations for the integration of the field 
lines are the same in each case, and form an array of regularly spaced grid points located near the lower boundary of the volume. 
Black field lines indicate (closed) lines that intersect the lower boundary twice, and red and green field lines represent field lines that 
leave the box through either the sides or top, with color indicative of polarity. 
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1 in Fig. [2]) 
is of lOkam 



a feature whic h was an- 
Qkamoto et alJ (l2QQ8k . Models 



the core volume shown 
ticipated by the analysis 
using the optimization method (McT, Wie, and Tha) contain 
between 4%-15% more energy than the potential solution, but 
possess more residual Lorentz forces than the current-field it- 
eration solutions. The magnetofrictional model (Val) has 4% 
more energy than the potential solution but has larger values 
ofV-B than the optimization or current-field iteration solu- 
tions. Based on the results summarized in Table [T] the excess 
magnetic energy (above the potential field) for this active re- 
gion could be anywhere from near zero to about 25% of the 
potential field energy. However, it is also possible that the ex- 
cess energy is significantly larger than 25% when taking into 
account the uncertainty associated with the inconsistency be- 
tween the boundary data and the force-free-model assumption 
(see S31). 

Because of these differences in the resulting NLFFF mod- 
els of AR 10953, we perform a goodness-of-fit test to de- 
termine which of the NLFFF models is the best approxima- 
tion to the observed coronal magnetic field. In the earlier 
study of lSchrijver et al.l ( 2008 ), we performed this test in both 
a qualitative and quantitative manner using EUV and X-ray 
imagery, provided respectively by the Transition Region and 
Coronal Explorer (TRACE) and Hinode/XRT instruments, by 
determining which model possessed field lines that were more 
closely aligned with the projected coronal loop structures vis- 
ible in the (two-dimensional) image plane. Models for which 
most field lines appeared to be aligned with loops were con- 
sidered good approximations to the actual coronal magnetic 
field. Locations where the field was noticeably sheared or 
twisted were of particular interest because such patterns are 
usually indicative of the presence of currents (which the mod- 
eling seeks to ascertain). More weight was typically given 
to regions connected to places at the photospheric boundary 
where Jz is found to be high, whereas coronal loops located 
in the periphery of the active region with footpoints located 
where Jz was lower were likely to be less sensitive to the 
presence of currents elsewhere in the active region. All such 
comparisons with coronal loops rest on the assumption that 
the plasma responsible for the emission is aligned with the 
coronal magnetic field and that this field is in a force-free 
state. 

For AR 10953, we overlaid field lines from all of the 
NLFFF models (as well as the potential field model) on top 
of the time-averaged Hinode/XRT image shown in Figure [T^, 
and used the same criteria listed above to qualitatively deter- 
mine the better-matching models. We subjectively judged the 
field lines in the Wh~ , Ami ~ , Am2~ , and Can" models to be 
more closely aligned with the XRT loops than any of the oth- 
ers. An overlay of field lines from the Wh~ model is shown 
in Figure[Tt. This judgement is based on good alignment with 
the tightly curved X-ray loops north of the sunspot (which is 
visible in the coaligned magnetogram of this region shown in 
Fig- [Hi), together with a reasonably good match of the loop 
arcade and fan structures to the south and west of the sunspot. 
This judgement is also based on side-by- side comparisons of 
field line overlays amongst the various candidate models (in- 



cluding the potential field model), from which a relative rank- 
ing was determined. The models listed above came out on top 
in both instances. 

With the aim of determining more quantitatively the best- 
fit model(s) for AR 10953, we also compared the model 
field lines to three-dimensional trajectories of loop paths. 
We are able to do this because AR 10953 was observed by 
the twin Solar Terrestrial Relations Observatory (STEREO) 
spacecraft, one of which leads the Earth in its orbit around 
the Sun, and the other of which trails the Earth. As part 
of the Sun Earth Connection Coronal and Heliospheric In - 
vestigation (SECCHI) instrument suite /Howa rd et aPbOOSh . 
each STEREO spacecraft contains an Extreme Ultraviolet 
Imager (EUVI). The angular separation of the two STEREO 
spacecraft at the time AR 10953 was on disk (of about 7°) 
was favorable for stereoscopically determining the three- 
dimensional trajectories of loops observed in the 17lA, 195 A, 
and 284A channels of EUVI. The coordinates of these loop 
trajectories were obtained by triangulating the positions of 
common features visible in p airs of concurrent EUVI im ages 
using the method described in lAschwanden et al] ( 2008k . 

Unfortunately, most of the loops visible in the three EUVI 
wavebands lie outside of the central region of AR 10953 
(Fig. [Tt), and thus do not overlap the region for which the 
vector magnetogram data are available (Figs.[T]l,e). The main 
reason is that loops located closer to the centers of the active 
regions tend to emit more in X-ray passbands than in EUV 
passbands. In addition, large loops at the periphery of active 
regions are generally easier to reconstruct with stereoscopy, 
while small loops in the centers of active regions are more 
difficult to discern from underlying bright features (such as 
moss) and thus cannot unambiguously be triangulated. How- 
ever, the outlying loops evident in AR 10953 should still sense 
the presence of currents in the center of the active region, 
due to Ampere's Law, and thus might be useful for quanti- 
tatively determining the best-matching NLFFF model for this 
active region. We infer that currents must be present in the 
AR 10953 corona for two reasons. First, most of the strong 
vertical currents in the Jz map are located in the central por- 
tion of the active region (as illustrated in Fig. [3]) and presum- 
ably flow upward into the corona. Second, field lines from the 
potential model do not qualitatively match the X-ray and EUV 
loops as well as field lines from the Wh~ , Ami ~ , Am2~ , and 
Can" models, which are our most nonpotential models and 
evidently contain currents strong enough to affect the trajec- 
tories of many field lines in the central portion of this active 
region (cf. Fig.O. 

To quantitatively compare the STEREO loops and the 
NLFFF-model field lines, we determine the (positive) angle (j) 
between the STEREO-loop and the model-field line trajecto- 
ries subtended at all STEREO-loop points lying inside the full 
320 X 320 X 256-pixel NLFFF computational domain. We then 
computed the mean of these angles, yielding for each model 
the domain-averaged misalignment angle metric {(j)) listed in 
Table [T] We find that, at least by this particular quantitative 
measure, none of the NLFFF models improve upon the value 
of {(j)) = 24° found for the potential field model, although 
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several models (including the qualitatively better- fitting mod- 
els discussed earlier) are comparable. We discuss reasons why 
none of the models improved upon the potential field metric 
for ((/)) in g21 



4. Discussion 



Given the boundary conditions produced using the data 
preparation process described in the various NLFFF algo- 
rithms converged to different solutions for the coronal field 
above AR 10953. A few of the models appear to match 
the loop structures in the Hinode/XRT image, but none of 
them were able to improve upon the potential field in their 
alignment with the three-dimensional loop trajectories in- 
ferred from STEREO/SECCHI-EUVI. In attempting to find 
a consensus model, we also applied the NLFFF algorithms to 
different boundary data generated using variants of the data 
preparation process. These variations, described in ^4.11 were 
run in parallel to those analyzed in ^ and also did not pro- 
duce a viable model. 

This inability to generate models that both qualitatively 
and quantitatively match the coronal loops paths is disap- 
pointing, especially given the generally successful applica- 
tion of these algorith ms to test cases with known solutions 
( Schrijver et al.l 120061) . including a solar-like test case with 
quasi-realistic forcing in the lower layers that was meant to 
approxi mate some of the fo rces acting in the solar chromo- 
sphere dMetcalf et al.l l2008l) . While we realistically expect 
the various methods to yield somewhat different solutions, we 
cannot fully ascribe the broad range of inconsistencies in the 
solutions solely to algorithmic differences. This causes us to 
examine the entire NLFFF modeling process from beginning 
to end, and in so doing we have identified several additional 
factors that likely also impact our ability to produce robust 
models. These factors are discussed further in ^4. 21 and ^4.31 

4.1. Data Preparation Variations 

We applied the NLFFF algorithms to boundary data pro- 
duced using eleven variations of the data preparation process, 
of which only one was outlined in Variations involved sub- 
stituting a different procedure to remove the 180° ambiguity 
of the measured transverse vector field, and/or using differ- 
ent versions of the standard preprocessing algorithm. In total, 
about 60 different NLFFF models for AR 10953 were cal- 
culated. (Not all algorithms were run on all of the available 
boundary data sets.) 

The first variant entailed using a different algorithm to re- 
move the 180° ambiguity inherent in the vector-magnetogram 
inversion process. Although there are in fact several algo- 
rithms to do this, we chose as an alternative to AZAM to 
employ t he automated Univer sity of Hawai'i Iterative Method 
(UHIM) dCanfield et al.lll993h because it has been used exten- 
sively in the literature and also sc ored highly amongs t other 
ambiguity resolution algorithms ( ^Metcalf et al.l |2006|) . We 
found that, while differences exist in, for example, field line 
trajectories near regions where the ambiguity was resolved 
differently, the volume-integrated metrics discussed in 0and 




I ' ' ' ' I ' ' ' ' I 
50 100 150 200 250 300 
pixels 



-10"^ -5x10^ 5x10^ 10"^ 
vertical current density [Mx cm"^ s"^] 



Fig. 3. — A map of the vertical component of the electric cur- 
rent density at the lower bounding surface as determined 
from Hinode/SOT-SP vector-field measurements (i.e., prior to 
preprocessing). The values of Bx and By, and therefore J^, 
outside of the region containing vector magnetogram data are 
unknown and have been zeroed out. Saturation of the color 
table is indicated by black or white hues. The pixel scale is 
580 km per pixel. 



shown in Table \T\ were largely similar for both the AZAM- 
and UHIM-disambiguated boundary data. 

The second variant involved a new version of the method 
used to preprocess the values of B\zq to make the boundary 
data more consistent with a force-free solution. Our standard 
scheme pivots and smooths the components of B\zq so that 
the integrated magnetic forces and torques in the overlying 
volume are reduced as much as possible, while also retaining 
some fidelity to the measured vector field. For AR 10953, we 
also experimented with a preprocessing scheme (described in 
Wiegelmann et al.ll2Q08|) that, in addition to the above, seeks 



to align the horizontal components of B\zq with fibrils seen in 
contemporaneous images of Ha. The motivation for this addi- 
tional preprocessing constraint is to produce boundary data as 
close as possible to the force-free field expected to exist at the 
chromospheric level (to which the Ha fibrils are assumed par- 
allel). We found, however, that using Ha-fibril information 
(observed by the Narrowband Filter Imager of Hinode/SOT) 
did not make a significant difference in the domain- averaged 
metrics used to characterize the various extrapolation models, 
although we intend to experiment further with this preprocess- 
ing scheme as it is somewhat new. 

The third var iant was to use thejnethod of preprocess- 
ing described in Fuhrmann et al. ( 2007), the goals of which 
are the same as the Wie gelmann et al. (200 6), but which uses 
a simulated annealing numerical scheme to find the optimal 
B\zQ field. As with the other variations, using this alternate 
preproc essing scheme did not m uch affect the resulting global 
metrics (iFuhrmann et al ] l2009h . 



4.2. Field-of-View Issues 

The Hinode/SOT- SP vector magnetogram data span only 
the central portion of the AR 10953, and thus do not cover all 
of the weaker field and plage that surround the active-region 
center. Here, as in the Schrijver et al. (2008) case, we chose to 
extend the NLFFF computational domain and embed the vec- 
tor data in a larger line-of-sight magnetogram. One benefit 
of such embedding is that it places the side and top bounding 
surfaces farther away from the center of the active region, in 
locations where the coronal magnetic field is presumed more 
potential and thus more consistent with the boundary con- 
ditions applied there. Another reason is that in earlier test 
cases using boundary data with known solutions (described 
in Schrijver et al. 2006), we found that enlarging the NLFFF 
computational domain improved the solution field in the cen- 
tral region of interest. We attributed this behavior primarily 
to the sensitivity of the final solution to the specified bound- 
ary conditions, and concluded that moving the side and top 
boundaries farther away from the region of interest improved 
the resulting models. 

However, there is an important difference between these 
earlier tests and the current case of AR 10953. In the 
ISchrijver et al.l (l2006f) study, vector data for the entire (en- 
larged) lower boundary were available, and thus the locations 
of currents penetrating the entire lower bounding surface, 
over both polarities, were known. In contrast, for AR 10953 
we have no information about currents located exterior to the 



region containing the Hinode vector magnetogram data, as 
shown in Figure [51 and consequently (as stated earlier) the 
horizontal components of ^l^o were set to zero in the region 
outside of the area containing Hinode/SOT- SP vector data. 
This is obviously not correct, but lacking any knowledge of 
actual horizontal fields there, this approach was presumed 
to be the least damaging. However, the lack of satisfactory 
results suggests that the decision to embed may not be as 
harmless as originally believed. 

The ability of the various NLFFF algorithms to find a valid 
solution ultimately depends upon how they deal with the cur- 
rents passing through the bounding surfaces of the computa- 
tional domain. Figure |4] shows maps of the current density 
integrated vertically through the models. It is evident from 
these images that algorithms based on similar methods result 
in models that look similar to each other, but also that there 
are stark differences between the locations of the strong cur- 
rents amongst the different classes of methods. 

It is interesting to note that for AR 10953, as for the 



Schrijver et al.l (120081) case, the solutions bearing the best re- 



semblance to the Hinode/XRT loops, and here were among 
the best at matching the STEREO-loop trajectories, were cal- 
culated using the current-field iteration method. This method 
differs from the others in that it uses values of Jz and a only 
in one of the polarities (the well-observed leading polarity, 
in the case of the best-fit models) from the lower boundary, 
while ignoring such measurements in the opposite polarity. In 
contrast, the optimization and magnetofrictional methods re- 
quire that information about currents be available across both 
polarities. 

We suspect that the Wheatland current-field iteration al- 
gorithm benefits from the additional space in the solution 
domain because fewer current-carrying field lines intersect 
the side boundaries (which causes their values of a to be 
set to zero). However, the Wiegelmann optimization algo- 
rithm, and the Valori magnetofrictional algorithm in particu- 
lar, perform better when applied to smaller volumes or when 
the weighting given to the peripheral boundary information 
is less than that applied to the Hinode vector magnetogram 
data. The bottom row of images in Figure H] shows that the 
Valori magnetofrictional algorithm has markedly different be- 
havior depending on the weighting of the peripheral boundary 
data. The differences are most striking in the area exterior to 
where the vector-magnetogram data is located. Restricting 
the computational domain to contain only the region over- 
lying the Hinode/SOT- SP field of view produces a solution 
with more intense currents and having fewer Lorentz forces 
((CWsin(9) = 0.19) and greater energy (E/Epot = 1.12) 
than the Val solution. Many of these problems caused by the 
embedding process are alleviated when vector magnetogram 
data are provided over a field of view that covers the locations 
of all relevant currents associated with the region of interest. 
For active-region studies, this often means capturing much of 
the trailing polarity, which is often more diffuse and extended 
than the leading polarity. 

We therefore conclude that vector magnetogram data of ac- 
tive regions for use by NLFFF modeling efforts need to span 
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Fig. 4. — Images showing the magnitude of the current density \J\ after integrating vertically through the computational domain for 
most of the models presented in Fig. [2l Algorithms using the same class of method tend to produce similar patterns, as evident in 
the top row (showing models produced using optimization algorithms) and in the middle two rows (showing models produced using 
Grad-Rubin- style current-field iteration algorithms). The bottom row illustrates three different versions of the Valori magnetofric- 
tional model, illustrating some of the effects associated with the process of embedding vector magnetogram data into line-of-sight 
magnetogram data to produce lower-boundary data. Shown are the Val model of Fig. [2] which weights more heavily the boundary 
data inside the Hinode/SOT-SP field of view, a model for which the lower-boundary data was weighted uniformly, and a smaller- 
domain model encapsulating only the volume overlying the Hinode/SOT-SP field of view. The integrated current map from the 
Am2~ model in Fig.[2lis almost identical to that of the Ami" model, and is not shown. 
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much of the area above which currents flow. Coverage of the 
more diffuse, traiUng-polarity fields is Hkely to be especiaUy 
important because of the tendency for the traihng-polarity 
field to contain the endpoints of many field lines that carry 
significant currents (due to the existence of such currents in 
the leading polarity, coupled with the assumption that many 
field lines connect the leading and trailing polarities within 
the active region of interest). 

On a related topic, we suspect that the STEREO-loop com- 
parison process described in 0is affected both by the prox- 
imity of the STEREO loops to the sidewalls of the NLFFF 
computational domain (where potential-field boundary condi- 
tions were applied) and by their lying outside of the region for 
which we have vector magnetogram data (Figs.[T]l,e). Conse- 
quently, one might not be surprised that the potential model 
bested the others in matching the STEREO loops, but the siz- 
able misalignment angle {(j)) of 24° for the potential model 
seems to suggest that even these outlying STEREO loops do 
carry some currents. 

In light of these issues, rather than using the STEREO- 
loop comparison as a discriminator between the collection of 
NLFFF models, we instead view the collectively poor mis- 
alignment angles by the NLFFF models as another indication 
that the region over which vector magnetogram data are avail- 
able needs to be enlarged. Although it is possible to enlarge 
the NLFFF computational domain (beyond what we have al- 
ready done) in order to include even more loops observed by 
STEREO, we again emphasize that the added benefit of do- 
ing so without additional vector magnetogram data would be 
minimal because of the lack of further information about cur- 
rents flowing through the lower boundary. Indeed, we applied 
the same current-field iteration method used for Wh~ to larger 
(512x51 2-pixel) boundary data produced using the same pro- 
cess described in ^ and found that the value of ((/)) for the 
identical volume used to compute the values of {(j)) in Table [T] 
remained unchanged. 

Lastly, we recognize that, when compared with stronger- 
field regions, the transverse field components Bt are not mea- 
surable with the same degree of certainty in weaker-field re- 
gions such as those likely to lie within the enlarged fields of 
view for which we are advocating. The findings presented 
here, however, suggest that the NLFFF modeling algorithms 
would benefit by having these vector magnetic field data avail- 
able, even if such data possess higher measurement uncertain- 
ties than the stronger fields found closer to the centers of most 
active regions. 

4.3. Boundary Data Inconsistencies 

In we described several conditions that the boundary 
data B \ s must satisfy in order to be consistent with a force- 
free magnetic field. However, these conditions are never guar- 
anteed to be satisfied on the full bounding surface S, which 
here consists of the vector and line-of- sight magnetogram 
data for the lower boundary combined with the potential field 
boundary conditions used for the remainder of the enclosing 
surface. To partially rectify this problem, we apply prepro- 
cessing to these data to thereby adjust the various compo- 
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Fig. 5. — (a) Scatter diagram illustrating the mismatch 
between the values of a implied by the preprocessed 
Hinode/SOT-SP boundary data B\zo for all points having 
Bz > 0, and the values of a for field lines in the Wh~ model 
intersecting these same points. For a consistent boundary con- 
dition where the a-correspondence relation is satisfied, the 
values of a on each field line in the Wh~ solution (which are 
taken from the negative-polarity end of the field line) would 
match the measured value of a found at the positive-polarity 
end. (b) The differential change d^/da in net flux ^ inte- 
grated over all points having a values greater than the ab- 
scissa. The unpreprocessed (red) and preprocessed (black) 
boundary data for AR 10953 are both shown. Data for which 
the ^-correspondence relation holds have d^/da = 0. 
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Fig. 6. — Censorship map for the Wh~ model. Pixels are gray 
in the positive polarity of Bz and are either black or white in 
the negative polarity. White pixels indicate the locations of 
field lines for which a was censored (set to zero); black pixels 
indicate locations where a 7^ 0. Contours of Bz are overlaid, 
with green and red contours located in positive and negative 
polarity regions, respectively. Contour levels are at ±{100, 
200, 500, 1000, 2000} Mx cm-^. The pixel scale is 580 km 
per pixel. 



nents of B\zq on S such that the boundary data are made more 
compatible with the equations the NLFFF algorithms seek to 
solve. 

Even after preprocessing, however, the boundary data can 
be shown to be incompatible with a force-free field. The Wh~ 
model, which is one of several models judged to match best 
on a qualitative basis, only uses the a values located in the 
negative polarity of the active region. However, the algorithm 
converged to a solution for which the corresponding a val- 
ues in the positive polarity do not match those indicated by 
the Hinode/SOT-SP data. Figure [5^ illustrates this problem. 
There, the a values in the Wh~ model from field lines that 
intersect the lower boundary in the positive polarity are plot- 
ted versus the a values at the same boundary points deduced 
from the preprocessed Hinode data. For consistent boundary 
data, these would be equal. The scatter evident in the figure 
indicates that the Hinode boundary data, even after prepro- 
cessing, are inconsistent with a force-free field. Additionally, 
the difference in the location of currents in the Wh~ and Wh+ 
models (and similarly in the Reg~ and Reg+ models), as evi- 
dent in Figure m may also indicate that the boundary data are 
inconsistent with a force-free solution. 

Figure [5j) illustrates this effect in a different way. This 
incompatibility can be illustrated by computing 



^{a) = / H{a - a) Bzdxdy, 
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(1) 



where H is the Heaviside step function, and Bz{x^y) and 
a'{x^ y) are, respectively, the flux density and value of Jz/Bz 
at each point on the preprocessed Hinode boundary map. The 
function ^{a) signifies the net flux in that subarea of the 
boundary map for which a' is larger than a certain threshold 
a. When the a-correspondence relation holds, the function 
^(a) thus possesses a derivative of zero because such cor- 
respondence requires, for any interval da, an equal amount 
of positive and negative flux passing through that subarea of 
the boundary map having values of a' between a and a^da. 
However, Figure[5]3 shows that d^ / da is nonzero over most a 
values for the preprocessed data used here, especially within 
the range —0.2 < a < 0.2 which corresponds to the a values 
possessed by about 80% of the area of the boundary map. For 
comparison, the figure includes the function d^ / da for the 
unpreprocessed dataset. 

The various methods deal with the lack of a correspon- 
dence in the boundary data in different ways. Current-field 
iteration methods allow the a-correspondence condition to be 
met by ignoring the values of a in one polarity. However, only 
limited uniqueness results have been found for this approach, 
and even existence r esults are limited to the case of an un- 
bounded domain f see lAmari et al ]l2006h . It is well known that 
the current-field iteration method fails to converge in some 
cases, and this may be due to the absence of a solution, or the 
absence of a unique solution. In Wheatland's implementation 
of this method, if the solution does not converge, values of a 
are censored (set to zero) in the polarity defining the currents 
going into the corona. The censorship is imposed at boundary 
points with | Bz \ less than a threshold value, and that value is 



increased as required. Additional censorship is also imposed 
such that field lines intersecting the side and top boundaries 
carry no current. In practice it is found that such reduction of 
the currents flowing into the domain can lead to convergence. 
The Wh~ model, for example, censored almost half of the 
values of a in the negative polarity (corresponding to 43% of 
the negative-polarity flux) before convergence was achieved, 
as illustrated in Figure [6l Valori's magnetofrictional method 
is prevented from relaxing past an equilibrium state in which 
the continual injection of inconsistencies into the model (at 
the boundaries) is balanced by their removal via diffusion. 
Wiegelmann's optimization method does not reach as well- 
relaxed of a force-free state as some of the other models, even 
though it disregards some of the boundary mismatches via the 
tapered nature of the weighting functions towards the edges 
of the model volume. 

There are several reasons why the boundary conditions 
used for this study (and other active region studies) might not 
satisfy the force-free consistency relations. The most conspic- 
uous reason is that the photospheric layers of the Sun, from 
which originate the Hinode/SOT-SP magnetogram data used 
here, do contain Lorentz, buoyancy, and pre ssure gradient 
force s and thus ar e not force-free to begin with ( Metcalf et al.l 
19951: lGarvl[2QQll) . Additionally, measurement uncertainties 
in the components of B\zq preclude accurate determinations 
of Jz (and thus a) on the lower boundary because of the need 
to take derivatives of the horizontal components of B\zq. An- 
other reason is that measurements of the current density nor- 
mal to the enclosing surface are unavailable over much of S 
due to the lack of vector magnetogram data above the photo- 
sphere. Another is that the modeling implicitly assumes that 
the boundary data span a planar surface, and do not take into 
account effects present in vector magnetograms such as the 
Wilson depression in sunspots and the broad range of line- 
formation heights across the line. Yet another is that the 
inversion techniques that produce the vector magnetogram 
measurements do not fully take into account the multiple 
components of thin, narrow strands of inte rleaved magnetic 
fields that chara cterize sunspo t penumbrae ([Title et al.lll993l: 
Beflot Rubio et al. 2004; ShimizuetaPllQQSb. We thus con- 
clude that the NLFFF modeling process needs to account for 
these intrinsic uncertainties in the boundary data, which in- 
clude everything from measurement uncertainties to the lack 
of knowledge about how to infer the magnetic field in the 
force-free region at the base of the corona from the observed 
photospheric field maps. 

5. Conclusions 

We have attempted to model the coronal magnetic field 
overlying AR 10953 by applying a suite of NLFFF algorithms 
to the photospheric vector field measured using Hinode/SOT- 
SR These data were remapped, embedded, and preprocessed 
in various ways in order to produce boundary data for this 
active region that were also consistent with the force-free 
assumption. From these boundary data, about 60 different 
NLFFF models were constructed. 

The resulting variations in these models prompted us to 



validate the results against images of coronal loops evident in 
EUV or X-ray images. The goodness of fit was first deter- 
mined in a qualitative manner by overlaying NLFFF-model 
field lines on Hinode/XRT imagery. This comparison indi- 
cated that some models contain field lines that are aligned 
with the observed loop structures. However, conclusive de- 
terminations of best-matching models, based solely on such 
overlays, remained difficult because of the indistinct nature of 
many coronal loops, especially those located near the center 
of AR 10953 where many of the currents are presumed to lie. 

We then turned to stereoscopic determinations of three- 
dimensional loop paths as a way to quantitatively assess 
the goodness of fit. This comparison was also incon- 
clusive, because the loops traced stereoscopically in the 
STEREO/SECCHI-EUVI observations were restricted to the 
outermost domain of the active region. This meant that those 
loops that did fall in the NLFFF computational domain lay 
close to the edge of the computational volume, where model 
field lines either leave the domain or run close to the side 
boundaries. We suspect this quantitative comparison was at 
least partially compromised by these effects, due to the model 
fields being sensitive to the way in which the side boundary 
information is incorporated and to their being located above 
the portion of the lower boundary for which Hinode/SOT-SP 
vector magnetogram data were not available. 

As exemplified by the qualitative and quantitative compar- 
isons presented here, we find that it remains difficult to con- 
struct and validate coronal magnetic field models of solar ac- 
tive regions that can reliably be used for detailed analyses of a 
quantitative nature. Our experience with modeling test cases 
with known solutions had shown that the various algorithms 
do work when given consistent boundary conditions. This led 
us to examine thoroughly the entire NLFFF modeling frame- 
work in order to identify problematic issues that impact our 
ability to build useful models of the solar coronal field. The 
results of this examination leave us with several possibilities. 
First, it may be that useful NLFFF extrapolations based on 
currently available signal-to-noise levels, preprocessing pro- 
cedures, fields of view, and observable fields are intrinsically 
infeasible. A second (and more hopeful) possibility is that 
NLFFF extrapolations need both much larger fields of view 
to better constrain the long field lines high over a region or 
to distant neighboring regions, and enough spatial resolution 
to resolve the spatial distribution of current densities on the 
boundaries. Third, NLFFF algorithms need to accommodate 
the fact that the boundary conditions contain (sometimes sig- 
nificant) uncertainties, either from the measurement process 
(e.g., signal-to-noise issues or inadequate resolution of the 
180° ambiguity), or from physical origins (e.g., variations in 
the line-formation height, or most prominently the non-force- 
free nature of photospheric vector magnetograms). 

The second possibility can be tested empirically. One way 
to do this with current codes and instrumentation is to obtain 
vector magnetic observations of a substantially smaller active 
region and its wide surroundings. This will place the side 
boundaries relatively farther away from the region of interest, 
while remaining compatible with the range and resolution of. 



13 



e.g., the Hinode/SOT-SP and with the Cartesian nature of the 
available modeling codes. 

To address the third possibility, we have several avenues 
available. Simple ways to account for boundary data un- 
certainties include introducing a position-dependent weight- 
ing function used in relaxation methods, or modifying the 
selection criteria for the a field in the current- field itera- 
tive method. Additionally, the preprocessing of the raw 
vector data needs to better approximate the physics of the 
photosphere-to-chromosphere interface in order to transform 
the observed photospheric field to a realistic approximation 
of the overlying near-force-free field at the base of the corona. 
One way to do that without resorting to more computationally 
intensive MHD models is to use the magne tohydro static con- 
cept (e.g., 'Wiegelm ann & NeukirchI 12006*) and approximate 
the stratifications for the flux tubes and their surroundings (or 
the strongly and weakly magnetic regions) separately. 

Finally, in light of our findings in this study and in con- 
sideration of the aforementioned goal of constructing models 
that provide useful estimates of physical quantities of inter- 
est, we thus recommend that a particular force-free extrapola- 
tion should not be considered a consistent model of an active- 
region corona unless the following indicators (at a minimum) 
are satisfied: (1) good alignment of modeled field lines to the 
coronal loops observed on the solar disk; (2) acceptable agree- 
ment of the a-correspondence relation by having similar val- 
ues of a at both ends of all closed field lines, and acceptable 
agreement with the boundary values of a from the data; while 
(3) still realizing low values of the NLFFF metrics (CW sin 6) 
and(|/i|). 
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